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Abstract. Over the coming decade, the observational samples available for studies of 
cluster abundance evolution will increase from tens to hundreds, or possibly to thou- 
sands, of clusters. Here we assess the power of future surveys to determine cosmological 
parameters. We quantify the statistical differences among cosmologies, including the 
effects of the cosmic equation of state parameter w, in mock cluster catalogs simulating 
a 12 deg^ Sunyaev-Zel'dovich Effect (SZE) survey and a deep 10^ deg^ X-ray survey. 
The constraints from clusters are complementary to those from studies of high-redshift 
Supernovae (SNe), CMB anisotropics, or counts of high-redshift galaxies. Our results 
indicate that a statistical uncertainty of a few percent on both flm and w can be 
reached when cluster surveys are used in combination with any of these other datasets. 



INTRODUCTION 

Because of their relative simplicity, galaxy clusters provide a uniquely useful 
probe of the fundamental cosmological parameters. The formation of the large- 
scale dark matter potential wells of clusters is likely independent of complex gas 
dynamical processes, star formation, and feedback, and involve only gravitational 
physics. The observed abundance of nearby clusters implies robust constraints 
on the amplitude erg of the power spectrum on cluster scales to an accuracy of 
~ 25% [1,2]. In addition, the redshift-evolution of the observed cluster abundance 
constrains the matter density Qq [3-5]. 

In order to be useful for these cosmological studies, the masses of galaxy clusters 
have to be known. Existing studies utilized the presently available tens of clusters 
with mass estimates [6,7], and as a result, were limited in their scope. The present 
samples, however, will likely soon be replaced by catalogs of thousands of inter- 
mediate redshift and hundreds of high redshift {z > 1) clusters. At a minimum, 
the analysis of the European Space Agency X-ray Multi-mirror Mission (XMM) 
archive for serendipitously detected clusters will yield hundreds, and perhaps thou- 
sands of new clusters with temperature measurements [8]. Dedicated X-ray and 
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SZE surveys could likely surpass the XMM sample in areal coverage, number of 
detected clusters or redshift depth. 

The imminent improvement of distant cluster data motivates us to estimate 
the cosmological power of future surveys. In particular, wc study the constraints 
provided by a 12 deg^ SZE survey [9], or by a deep 10'^ deg^ X-ray survey. Our 
primary goals are (1) to quantify the accuracy to which various cosmological models 
can be distinguished from a standard A Cold Dark Matter (ACDM) cosmology; and 
(2) to contrast constraints from clusters to those available from CMB anisotropy 
measurements, from high-redshift SNe, or from high-redshift galaxy counts [10-12]. 

The galaxy cluster abundance provides a natural test of models that include a 
dark energy component with an equation of state parameter w = p/p ^ —1 [13-15]. 
The value of w directly affects the linear growth of fluctuations, and the angular 
diameter distance (and hence the SZE decrement and the X-ray luminosity) to 
individual clusters. We restrict our analysis to a flat universe, and focus on the 
following four parameters: the matter density Q^; the equation of state parameter 
w (assumed to be constant); the Hubble constant h = i^o/lOO km s"^ Mpc~^' and 
the amplitude of the power spectrum on 8/i"^Mpc scales, cTg. A broader range of 
parameters, including open/closed universes, and evolving w{z), will be examined 
in future work [16]. Details of the study described here can be found in [17]. 

MODELING DETAILS 

General Approach 

To quantify the power of a future cluster survey to distinguish cosmologies, we 
utilize the following approach: 

1. We pick a fiducial cosmological (ACDM) model, with the parameters 
{Q\,Q„^,h,(Ts,n) = (0.7,0.3,0.65,0.9,1), based on present large scale struc- 
ture data [18]. We assume this model describes the "real" universe. 

2. In the fiducial model, we compute the abundance of clusters dNf[^/dz as a 
function of redshift in a specific (SZE or X-ray) survey. This simulates the 
dataset that will be available in the future for cosmological tests. 

3. We vary parameters of our model, and recompute the cluster abundance 
dN^est/dz as a function of redshift in this new "test" cosmology. In each 
test cosmology, we set the value of erg by requiring the local cluster abundance 
at redshift z to match the value in the fiducial cosmology. 

4. We compute the likelihood of observing the redshift distribution dNfi^/dz if the 
true distribution were dNtest/dz. We utilize both the normalization and shape 
of the distributions, by combining the standard Poisson and Kolmogorov- 
Smirnov tests. 

5. We repeat steps 3 and 4 for a wide range of values of w, flm, and h. 



Predicting Cluster Abundance and Evolution 



The fundamental ingredient of this approach is the cluster abundance, given a 
cosmology and the parameters of a survey. In this study, we utilize the "universal" 
halo mass function found in a series of recent large-scale cosmological simulations 
[19]. Following these simulations, we assume that the comoving number density of 
clusters at redshift z with mass M is given by 
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where ctm is the present day r.m.s. density fluctuation on mass-scale M [20], 
is the linear growth function, and po is the present-day mass density. The directly 
observable quantity is the number of clusters with mass above Mmin at redshift 
z±dz/2 in a solid angle dfi: 
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where dV/dzdfl is the cosmological volume element, and Mmin(^) is the limiting 
mass of the survey, as discussed below. Equations 1 and 2 reveal that the cluster 
abundance depends on cosmology through several quantities: (1) the growth func- 
tion Dz] (2) the volume element dV/dzdfl; (3) the power spectrum ctm; (4) the 
mass density Pq; and (5) the survey mass threshold Mmin. The first four of these 
dependencies are well-determined, once the parameters of the cosmological model 
and the power spectrum are specified (in particular, note that the abundance is 
exponentially sensitive to the growth function). The scaling of the limiting mass 
with cosmology depends on the specific survey. 



Cluster Surveys 

In this study, we examine two specific surveys, in which clusters are detected 
through either their SZ decrements, or X-ray fiuxes. In practice, the only survey 
details we utilize in our analyses are the virial mass of the least massive, detectable 
cluster (as a function of redshift and cosmological parameters), and the solid angle 
of the survey. The SZE survey we consider is that proposed by J. Carlstrom and 
collaborators [9] . This interferometric survey will detect clusters more massive than 
~ 2 X 10^^ Mq, nearly independent of their redshift, and will cover an area of 12 dcg^ 
in a year. The relatively small solid angle of the survey will allow cluster redshifts 
to be determined by deep optical and near infrared foUowup observations. The X- 
ray survey we consider is similar to a proposed Small Explorer mission, called the 
Cosmology Explorer, spear-headed by G. Ricker and D. Lamb. The survey depth 
is 3.6 X 10^ cm^s at 1.5 keV, and the coverage is 10^ deg^ (approximately half the 
available unobscured sky). We focus on clusters which produce 500 detected source 



counts in the 0.5:6.0 keV band, sufficient to reliably estimate the emission weighted 
mean temperature. The X-ray survey could be combined with the Sloan Digital 
Sky Survey (SDSS) to obtain redshifts for the clusters. 

The most important aspect of both surveys is the limiting halo mass 
M-^ia{z,flrn,u!,h), and its dependence on redshift and cosmological parameters. 
For an interferometric SZE survey, the relevant observable is the cluster visibility 
V, which is proportional to the total SZE flux decrement. The detection limit as 
a function of redshift and cosmology for this survey has been studied using mock 
observations of simulated galaxy clusters [21]. In the X-ray survey, M^ain follows 
from the cluster X-ray luminosity - virial mass relation [22]. Illustrative examples 
of the mass limits in both surveys are shown in Figure 1, both for ACDM and 
for a. w = —0.5 universe. The SZE mass limit is nearly independent of redshift, 
and changes little with cosmology. As a result, the cluster sample can extend to 
z ^ 3. In comparison, the X-ray mass limit is a stronger function of w, and it rises 
rapidly with redshift. For the X-ray survey considered here the number of detected 
clusters beyond z 1 is negligible. The total number of clusters in the SZE survey 
is ~ 200, while in the X-ray survey, it is ~ 2, 000. 
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FIGURE 1. Limiting cluster virial masses for detection in a 10^ deg^ X-ray survey (upper pair 
of curves) and in a 12 deg^ SZE survey (lower pair of curves). The solid curves show the mass 
limit in our fiducial flat ACDM model, with w = —1, il,n — 0.3, and h — 0.65, and the dotted 
curves show the masses in the same model except with w = —0.5. 



RESULTS AND CONCLUSIONS 



Our results are summarized by the likelihood contours in the — w plane, 
shown in Figures 2 and 3 for the SZE and X-ray surveys, respectively. Figure 2 
shows three different cross-sections of constant total probability in the SZE survey, 
at fixed values of h (0.55,0.65, and 0.80) in the investigated 3-dimensional flm, h 
parameter space. These diagrams demonstrate that the constraints on Qm are at the 
~ 10% level, while w remains essentially unconstrained {w ^ —0.2). Nevertheless, 
the narrowness of the contours in Figure 2 implies that the SZE survey can yield 
accurate constraints on w if combined with other data. We find that the differences 
among cosmologies in the SZE survey are driven nearly entirely by the growth 
function Dz- This results from the cluster sample extending to high redshifts {z > 
1), where the growth functions in different cosmologies diverge rapidly. This makes 
the SZE sample especially useful. For comparison, galaxy counts at 2; ~ 1 probe 
mostly the cosmological volume, and constitute and independent test from clusters 
[12]. Note that our constraints scale weakly with h: this arises from the weak 
dependence of the power spectrum on h. 

Figure 3 shows constraints in the X-ray survey for a fixed h = 0.65. The increased 
number of clusters translates to significantly narrower contours compared to the 
SZE survey. The orientation of the contours remains similar, but we find that the 
cosmological sensitivity arises nearly entirely from the mass limit Mmm- Indeed, the 
X-ray flux is more sensitive to cosmology than the SZE decrement (cf. Fig. 1), and 
the X-ray sample extends only out to 2; ~ 1, where the growth functions are less 




FIGURE 2. Contours of 1, 2, and 3(7 likelihood for different models when they are compared to 
a fiducial flat ACDM model with fi^ ~ 0.3 and h = 0.65, using the SZE survey. The three panels 
show three different cross-sections of constant total probability at fixed values of h (0.55,0.65, 
and 0.80) in the investigated 3-dimensional flm, w, h parameter space. 



divergent. Also shown in Figure 3 are constraints expected from CMB anisotropics 
and from high-z SNe. The dashed curves correspond to the CMB constraints (a 
±1% determination of the position of the first Doppler peak); the dotted curves 
to the constraints from SNe (a ±1% determination of the luminosity distance to 
z = 1). As these curves show, the constraints from the CMB and SNe data are 
complementary to the direction of the parameter degeneracy in cluster abundance 
studies, making the cluster surveys especially valuable. 

Our findings suggest that cluster surveys will lead to tight constraints on a com- 
bination of Qrri and w, especially valuable because of their high complementarity to 
constraints from CMB anisotropics, magnitudes of high-z SNe, or counts of high-z 
galaxies. In combination with either of these data, clusters can determine both Qrn 
and w to a few percent accuracy. We have focused primarily on the statistics of 
cluster surveys: further work is needed to clarify the role of systematic uncertain- 
ties, arising from the cosmology-scaling of the mass limits and the the cluster mass 
function, as well as our neglect of issues such as galaxy formation in the lowest 
mass clusters. 

We thank J. Carlstrom and the COSMEX team for providing access to proposed 
instrument characteristics. ZH acknowledges support from a Hubble Fellowship. 




FIGURE 3. Likelihood contours for a fixed h — 0.65 in the X~ray survey. Also shown are 
combinations of w and Qm that keep the spherical harmonic index £ of the first Doppler peak in 
the CMB anisotropy data constant to within ±1% (dashed lines); and combinations that keep 
the luminosity distance to redshift z — I constant to the same accuracy. 
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